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1. Introduction 

Understanding the deconfinement phase diagram of QCD is one of the major challenges of 
particle physics. In this talk I will report on a description of the deconfinement phase transition 
at zero baryon density within the Hamiltonian approach to Yang-Mills theory in Coulomb gauge, 
extending the previously developed variational approach for the Yang-Mills vacuum lHJ to finite 
temperatures. This approach has been quite successful in understanding the infrared properties of 
Yang-Mills theory: One finds a gluon energy which is infrared divergent [|]] and a linearly rising 
quark potential both being signals of confinement. In addition, one finds an infrared enhanced 
running coupling constant with no Landau pole [||], a perimeter law for the 't Hooft loop [Q] and, 
within an approximate Dyson-Schwinger equation, an area law for the spatial Wilson loop [g]. 
Also the topological susceptibility was found in accord with the lattice data [^]. Given the success 
of this approach in the vacuum sector we can also expect a decent description of the deconfinement 
phase transition. The outline of my talk is as follows: I will first review the basic ingredients of 
the Hamiltonian approach to Yang-Mills theory in Coulomb gauge. Within this approach I will 
then study the grand canonical ensemble in a mean-field type approximation and investigate the 
deconfinement phase transition. 

2. Hamilton approach to Yang-Mills theory 

The Hamilton approach to Yang-Mills theory is based on the canonical quantization in Weyl 
gauge Ao = 0, which yields the Hamiltonian 

H= l - j<£x(n 2 (x)+B 2 (x)) , (2.1) 

where IT(x) = — i8 / 8A(x) is the canonical momentum operator. Due to the use of the Weyl gauge 
Gauss' law escapes the quantum equations of motion and has to be imposed as a constraint to the 
wave functional 

DU\i/[A}=0. (2.2) 

Here D denotes the covariant derivative in the adjoint representation of the gauge group. The op- 
erator in Gauss' law DTI is nothing but the generator of (time-independent but space-dependent) 
gauge transformations and Gauss' law ensures that (in the absence of matter fields) the wave func- 
tional must be gauge invariant. To respect gauge invariance one can work with explicitly gauge 
invariant wave functionals, which has been pursued mainly in 2 + 1 -dimensions but which be- 
comes exceedingly involved in 3 + 1-dimensions. A more convenient way is to fix the gauge and 
resolve Gauss' law explicitly. For this purpose Coulomb gauge is convenient and after gauge fixing 
one finds the following Hamiltonian 

H=-J d 3 x(J A 1 U(x)J A -n(x)+B 2 (x))+H c , (2.3) 
where Ja = Det(— Dd) is the Faddeev-Popov determinant and 

Hc = f / d3xd3 y J A l P a (x)JAF ab (x,y)p b (y) (2.4) 
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Figure 1: Comparison of the results of the variational calculation and lattice data for the gluon propagator. 

is the so-called Coulomb Hamiltonian. Here, p a = —f abc j^ h Yl c is the color charge density of the 
gluons. The gauge fixed Hamiltonian is highly non-local due to the presence of the Faddeev-Popov 
determinant and the Coulomb term H c . Although the gauge fixing gives rise to a more complicated 
(non-local) Hamiltonian it has the advantage that the gauge invariance has been taken care of once 
and for all and that after Coulomb gauge fixing any wave functional depending on the transversal 
gauge field only lies in the physical Hilbert space. Fortunately it turns out that in the gluon sector, in 
particular for the investigation of the infrared properties, the Coulomb term can be ignored, which 
I will do in the following. 



With the Coulomb gauge fixed Hamiltonian (2.3) the Yang-Mills Schrodinger equation has 
been solved by a variational principle using Gaussian type ansatze for the vacuum wave functional. 
The approach developed in our group differs from previous attempts [Q, ^ by the ansatz for the 
vacuum wave functional, by the full inclusion of the Faddeev-Popov determinant and by the renor- 
malization. Our ansatz for the vacuum wave functional reads 

1 



W [A] 



.-1/2 



exp 



A co A 



(2.5) 



where CO is a variational kernel, which is determined by minimizing the energy (y\H\y)- The pre- 
exponential factor has the advantage that it cancels the Faddeev-Popov determinant in the scalar 
product of the Coulomb gauge fixed wave functionals. Furthermore, for this wave functional the 
gluon propagator is given by 

(AA) = (2co)- 1 , (2.6) 

which identifies the Fourier transform of CO as quasi-gluon energy. Figure [T] shows the result of 
the variational calculation for the gluon propagator in comparison with the lattice data. As one 
observes, the cross feature of the lattice data are well reproduced. In particular, in the deep infrared 
regime our variational calculations agree perfectly with the lattice data. There are deviations in 
the mid-momentum regime, which can be substantially reduced by using non-Gaussian wave func- 
tionals, Ref. [§]. What is also remarkable is that the lattice data can be nicely fitted by Gribov's 
formula 

co(p) = yJp 2 +M 4 /p 2 (2.7) 

with the mass M = 0.88 GeV. 
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3. The grand canonical ensemble 

Since the gluons have vanishing chemical potential their grand canonical ensemble is defined 
by the density operator 

D = exp (-/?#), (3.1) 



where H is the gauge fixed Hamiltonian (2.3) and j8 = X/ksT is the inverse temperature. By means 
of the density matrix thermal averages are defined by 

( \ - Tr ( D -'-) H2) 
i~')T—fjD-> M 

where the trace can be calculated, in principle, in any complete basis. However, we will necessarily 
have to introduce approximations and then the choice of the basis matters. We will choose an 
optimal basis by exploiting the variational principle. Inspired by the zero temperature calculations, 
see Eq. ( |2.5[ ), we will choose our basis in the form 

\~k)= M 1 (3-3) 

where the states \k) are defined by the Fock basis 

|0), J{p)\0), a\p { )a\p 2 )\0), ... (3.4) 
obtained by decomposing the gauge field in terms of creation-annihilation operators 

A{k) = ^±={a(k)+a\-k)) (3.5) 
V 2co ( k ) 

with a so far arbitrary kernel co(k) and the vacuum state, being defined by a(k)\0) = 0, is the 
Gaussian 

(A\0)=cx P (-^Jacoa] . (3.6) 



Inserting this state into Eq. (3J) we find that the state |0) is nothing but the trial ansatz ( [Q| ) for 
the zero temperature variational calculation. However, in the present case the kernel co(k) is not 
determined by minimizing the energy but at the moment an arbitrary kernel. 

We cannot treat the full density operator, Eq. (3. 1 ). Therefore we replace in the density opera- 
tor the full Hamiltonian ( |2.4| ) by a single-particle one 

D = exp(-ph) , h = j dpe{p)a\p)a{p). (3.7) 

As a consequence, for the thermal averages (3.2) Wick's theorem applies, which tremendously 
simplifies the calculations of thermal expectation values. With the density matrix (3.7) one finds 
the usual Bose occupation numbers 

{a\p)a{p)) = n{p) = (exp(j3e(p)) - l)- 1 , (3.8) 

while the finite-temperature gluon propagator is given by 

(AA) T = (l + 2n)/(2co). (3.9) 
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It differs from the zero-temperature propagator (1.6) only by the presence of the finite-temperature 
occupation numbers n(p). In the limit of a vanishing temperature T — > the n(p) vanish and 



Eq. ( p.9| ) reduces to the zero-temperature propagator, Eq. ( 



With the density operator (3.7) at hand, we can straightforwardly calculate the entropy 



-k B TvD\nD 



(3.10) 



and the free energy 



F[e,G)] = (H) T -TS. 



(3.11) 



A comment is here in order: We cannot calculate in here the free energy from the partition function 
Z = TrD since the density matrix D (3^7) is not yet known. Rather we calculate the free energy by 
taking the thermal expectation value of the full Hamiltonian. This has in addition the advantage that 
also two-body correlations are included, which cannot be captured by the single-particle density 
matrix (|3~7|). 



4. The deconfinement phase transition 



So far, we have two unknown kernels. The single-particle energies e(p) occurring in the 
single-particle density operator (3/7) and the kernel (o(p) occurring in the vacuum state ( ^^ ) of 



our Fock basis (3.4), (3.6). We now determine these kernels by minimizing the free energy ( 3.11 ). 
Since n(p) is a monotonous function of e(p) we can vary F[e,n] with respect to n(p), which yields 



the relation [11] 



e(p) = co(p). 



(4.1) 



This relation is not surprising given the form of the thermal gluon propagator (3.9), which identifies 
(o(p) as quasi-gluon energy. Let us stress, however, that if one goes beyond the present approxi- 



mation and includes also the Coulomb term ( |2.4| ), one obtains still a linear relation between e(p) 
and (o(p), however, the proportionality factor differs from one by an additional loop integral 
Finally, variation with respect to the choice of our basis, i.e. with respect to (o(p) yields the gap 
equation 



(0 1 {p)=p 2 + X 2 {p)+I[n\. 



(4.2) 



where %(p) denotes the ghost loop calculated from the finite-temperature gluon propagator and 
I[n] denotes the tadpole of the transversal spatial gluons, see Fig. ||. The gap equation ( p~2| ) has 
to be solved together with the Dyson-Schwinger equation for the ghost form factor d(p). These 
equations can be solved analytically in the infrared by power law ansatze 



A B 



(4-3) 



Assuming the horizon condition d 1 (p = 0) = one finds at zero temperature from the ghost DSE 
the sum rule 



a = 2/3 + 2 -d, 



(4-4) 
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Figure 2: Tadpole term I[n]. 
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Figure 3: Infrared exponent j3 of the ghost form 
factor as a function of temperature. 



where d is the number of spatial dimensions. With this sum rule one finds then from the gap 
equation the following solutions for the infrared exponent j8 



d 
d 



3 : 
2 : 



y3 = l.oo, £ 

j8 =0.50. 



0.795 



(4.5) 



The same exponents are extracted numerically from the self-consistent solution of the equations of 
motion. 

At arbitrarily finite temperature an infrared analysis is impossible due to the fact that the gluon 



energy co(p) occurs in the occupation number ( p3§) in exponential form. One can, however, carry 
out the infrared analysis at infinite temperature, where the Bose occupation numbers given by 
Eqs. (0), (go]) reduce to 

T^oo: n(p) ->■ (li(o(p)y 1 . (4.6) 



One finds then the same sum rule ( |4.4| ) as at zero temperature and a single solution with 

T^oo: d = 3: j3 = 0.5 , a = 0. 



(4.7) 



Figure |3| shows the infrared exponent of the ghost form factor j8 as function of temperature. At 
low temperature one finds the two solutions for the infrared exponents obtained from the zero- 
temperature infrared analysis. These exponents stay more or less constant as the temperature in- 
creases up to a critical temperature, where the two solutions for the infrared exponent j3 merge 
to a single one, which with increasing temperature approaches the value /3 = \ predicted by the 
infrared analysis in the high temperature limit. Figure || shows the numerical solution for the ghost 
form factor d(p) and the gluon energy co(p) for various temperatures. The obtained results are in 
agreement with the predictions of the infrared analysis. At zero temperature both d{p) and (o(p) 
are infrared diverging. Above T c the ghost form factor is still infrared diverging but its infrared ex- 
ponent is halved while the gluon energy becomes infrared finite in accordance with the prediction 
of the infrared analysis ([OJ). Furthermore above T c the plateau value of the gluon energy in the 
infrared increases with the temperature. This is seen in Fig. |5[ where we show the gluon energy at 
an infrared cut-off Xjr. co(Xm) = m(T) jumps at the deconfinement phase transition from a very 
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Figure 4: Zero- and finite-temperature solutions for the ghost form factor d(p) (left panel) and the gluon 
kernel co(p) (right panel). 
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Figure 5: The effective gluon mass co(X]r) as a Figure 6: Critical behavior of the effective gluon 
function of temperature. mass for T — > T c . 



large to a small value and afterwards increases linearly with the temperature. Zooming into the 
transition regime, which is done in Fig. ^ one can extract the critical exponent of co(Xir) = m(T) 
defined by 



in. 



(T)~(T/T C -K 



(4.8) 



and obtains the value K = 0.37. This value compares well with the value K = 0.41 obtained in 
Ref. [|12|] in a phenomenological quasi-particle model for the gluons using input data from the 
d = 3 Ising model, which is in the same universality class as SU(2) gauge theory. Fitting our scale 
at the value of the Gribov mass M = 0.88 GeV, see Eq. ([L7|), determined in Ref. Jl3| ] on the lattice, 
we find a critical temperature for the deconfinement phase transition of T c = 275 . . . 290 MeV. This 
is in reasonable agreement with the lattice value of T c = 295 MeV. This result is also in reasonable 
agreement with the value T c ~ 270 MeV obtained from the effective potential of the Polyakov loop 



calculated in Ref. [ |14| ] within the Hamiltonian approach. Let me also mention that recently the 
Hamilton approach to Yang-Mills theory in Coulomb gauge has been also used to study glueballs 



at finite temperature [15]. The results obtained in the gluon sector at finite temperature are quite 
encouraging for an extension of the Hamilton approach to full QCD at finite temperature and baryon 
density. 
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